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ABSTRACT 

The Chandra X-ray Observatory recently discovered an X-ray point source near the center of Cassiopeia 
A, the youngest known Galactic supernova remnant. We have conducted a sensitive search for radio 
pulsations from this source with the Very Large Array, taking advantage of the high angular resolution 
of the array to resolve out the emission from the remnant itself. No convincing signatures of a dispersed, 
periodic source or of isolated dispersed pulses were found, whether for an isolated or a binary source. 
We derive upper limits of 30 and 1.3 mJy at 327 and 1435 MHz for its phase-averaged pulsed flux 
density. The corresponding luminosity limits are lower than those for any known radio pulsar with age 
less than 10 4 years. Our search sensitivities to single pulses were 25 and 1.0 Jy at 327 and 1435 MHz. 
For comparison, the Crab pulsar emits roughly 80 pulses per minute with flux densities greater than 100 
Jy at 327 MHz and 8 pulses per minute with flux densities greater than 50 Jy at 1435 MHz. These limits 
suggest that Cas A belongs to the growing population of young neutron stars that are radio quiet 



1. INTRODUCTION 

Cassiopeia A, with an age of only 320 years, is the 
youngest known Galactic supernova remnant. The super- 
nova that gave rise to this remnant was probably discov- 
ered in 1680, when Flamsteed (1725) obs erved a 6th mag - 
nitude star within 13' of Cas A's location ( Ashworth 1980] ) . 
Because its abundances of heavy elements are consistent 
with those expected from the explosion of a massive star 
(Hughes ct al. 2000), Cas A is thought to be the remnant 
of a Type II supernova. Given the stellar initial mass func- 
tion, a large fraction of Type II supernovae are expected to 
leave behind neutron stars. While many searches at vari- 
ous wavelengths for the compact remnant in Cas A have 
been undertaken, only upper limits had been established 
until the recent Chandra discovery of an X-ray poin t source 
(hereafter X PS) near the center of the remnant (Lamm] 
baum 1999[). This XPS w as later confirmed in archival 
ROSAT (|Aschenbach 199S| ) and Einstein ( |Pavlov fc Za 



vlin 1999) images. The XPS is within 5" of the expansio n 
center of the remnant (van den Bergh fe Kampcr 1983|). 



For a remnant distance of 3.4+q 3 kpc, calcula ted through 



radia l velocity, proper motion, and age data (Reed et al 



1995), this offset corresponds to a maximum transverse ve- 
locity of roughly 200 km s~\ a t ypical velocity for young 



pulsars flCordes fc Chcrnoff 1998j ). 

Pavlov et al. (2000) find that the X -ray spectrum of the 
XPS is fit well by a power law with index 2.6 — 4.1, with 



luminosity (2 — 00) x 10 eigs s . Similarly, ChakiabaiLy 



I al. (2001) LsLiinaLL a puwti law indtx uf 2.0 3.0 and 

luminosity of (2 16) x 10 34 ergs s~L These indices arc 

much steeper than the X-ray spectra of Active Galactic 



Nuclei (AGN). Furthermore, given the measured surface 
density of AGN, the probability of finding an AGN that 
close to the center of the SNR is negligible. For these 
reasons, Pavlov et al. (2000) conclude that the XPS is 
almost certainly associated with the Cas A remnant. The 
fitted power law indices for the XPS are also much steeper 
than those measured for the 6 young X-ray pulsars hav- 
ing ages less tha n 10 4 years, whose spectra l indices range 
from 1.1 — 1.7 ( |Chakrabarty et al. 2001 ). In addition, 
the estimated XPS luminosities are only marginally consis- 
tent with measured young pulsar X-ray luminosities, which 
range from 3 x 10 35 to 4 x 10 36 ergs s . It is also signif- 
icant that, unlike 5 of the 6 young pulsars with ages less 
than 10 4 years, Cas A shows no evidence for a synchrotron 
nebula surrounding the point source. 

While there is much evidence for and against the neu- 
tron star nature of the XPS, the definitive proof would 
be the detection of pulsations. No X-ray pulsations have 
been found in data taken with the Chandra High Resolu- 
tion Camera. However, a time-tagging problem caused the 
timing accuracy to be degraded from 16 /is to w 4 ms, al- 
lowing searc hes for only broad pulses w ith periods greater 
than 20 ms (Chakrabarty et al. 2001). The detection of 
radio pulsations, however, would be just as conclusive. For 
this reason, we have conducted a search for radio pulsa- 
tions from the XPS using the Very Large Array (VLA) at 
frequencies of 327 and 1435 MHz. There have be en several 
previous searches for a radio pulsar in Cas A (|Davies fc 
Large 197C|, |Seiradakis fc Graham 1980| , |Woan fc Duffctt 

19981 ) 



Smith 1992, Lorimer et al. 



Our VLA search is 
many times more sensitive, as having an accurate position 
for the XPS allows us to resolve out the strong emission 
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from the remnant, which with 327- and 1435-MHz flux 
densities of 6400 and 2000 Jy is one of the strongest radio 
sources known. 

2. DATA 

We observed the XPS on three days (23 October 
1999, 31 October 1999, and 28 November 1999) with 
the phased-array VLA in B array, where the maximum 
baseline is 11.4 km. The J2000 coordinates of the XPS 
are RA=23 h 23 m 27 s 94, DEC=+58°48 / 42", with positional 
uncertainty of 1" (Tananbaum 1999). As the 327- and 
1435-MHz synthesized beamwidths are approximately 20 
and 4 arcseconds respectively, one beam at each frequency 
was sufficient for this search. During each session, the XPS 
was observed at center frequencies of 327 and 1435 MHz, 
resulting in three pointings of length 20, 40, and 50 min- 
utes at 327 MHz, and four pointings of length 40, 40, 45, 
and 50 minutes at 1435 MHz. During each observing ses- 
sion, we observed the known pulsar PSR B0355+54, which 
has a period of 154 ms and a dispersion measure (DM) of 
57 pc cm" 3 . A blank-sky pointing and the flux and phase 
calibrator J2350+646 (with flux densities of 27 and 5 Jy 
at 327 and 1435 MHz) were also observed. 

We used the VLA's High Time Resolution Processor 
(Moffett 1997) to record dual circular polarization data at 
a sampling frequency of 9.22 kHz. Each scan was started 
on a 10-second tick, tied to the VLA's hydrogen maser, 
compared to Universal Time through the GPS. At 327 
MHz, 14 x 0.125-MHz frequency channels covered a total 
bandwidth of 1.75 MHz. These channels were spaced non- 
contiguously to avoid known interference frequencies and 
covered a total bandwidth range of 3 MHz. At 1435 MHz, 
14 x 4-MHz channels covered a total contiguous bandwidth 
of 48 MHz. (As only 50 MHz of bandwidth is available, 
some channels were duplicated) . The gain in each of the 
28 channels (2 polarizations and 14 frequencies) was cal- 
culated from observations of the calibrator source. 

3. ANALYSIS 

For each of the seven XPS data sets, we first summed 



right md left circular polarizations and then dedispersed 



the data over a range of trial DMs. The Taylor and Cordes 
(1993) model predicts a maximum DM of approximately 
75 pc cm" 3 in the direction of the XPS. We dedispersed 
the data with 60 trial DMs from to 300 pc cm" 3 to allow 
for an underestimate or a contribution from Cas A. These 
60 DMs were spaced nonuniformly, with greater spacing 
at higher DMs, where the temporal smearing due to dis- 
persion across an individual frequency channel is larger. 
Each dedispersed time series was Fourier transformed with 
32 Mpoint and/or 16 Mpoint Fast Fourier Transforms 
(FFTs), depending on the length of the data set. Shorter 
FFTs of 8 and 4 Mpoints were done to increase our sensi- 
tivity to some binary pulsars. The resulting power spectra 
were harmonically summed (up to a maximum of 1 6 har- 
monics) to incr ease our sensitivity to narrow pulses (Lyne 
& Smith 1998). All features in the resulting power spec- 
tra with signal-to-noise above a threshold of seven were 
saved as pulsar candidates. For each candidate, the corre- 
sponding dedispersed time series was folded at the period 
corresponding to the fundamental harmonic. These pro- 
files were inspected and candidate lists from different days 



and at different frequencies were compared. While many 
periodic signals were detected, they were predominantly 
at low DM, had noiselike pulse profiles, and did not ap- 
pear on multiple epochs or at multiple frequencies, and 
were therefore attributed to radio frequency interference 
(RFI). Figure 1 shows the folded pulse profiles resulting 
from applying the search algorithm to PSR B0355+54. 
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Fig. 1.— Folded pulse profiles for PSR B0355+54, which has a 
period of 156 ms and DM of 57 pc cm~ 3 , at 327 MHz (left) and 
1435 MHz (right). Note that these profiles are for 5 and 2 minutes 
of data at 327 and 1435 MHz, whereas our longest XPS search scans 
were 50 minutes at these same frequencies. 



An independent Fourier analysis of the same data, car- 
ried out at the New Mexico Institute of Mining and Tech- 
nology, was also unsuccessful. This analysis was sensitive 
to periods >10 ms and covered a smaller DM range, but 
with a finer DM spacing. Three filter-bank channels with 
high levels of RFI were rejected, decreasing the number of 
spurious candidates significantly with only a 12% sensitiv- 
ity loss. Power spectra were harmonically summed up to a 
maximum of 32 harmonics. Final pulsar candidates were 
evaluated by comparing average profiles from two halves of 
the bandpass and two halves of profiles folded over a two 
period span. These were compared with the corresponding 
profiles for non-dedispersed data. 

The Crab pulsar, currently the youngest known radio 
pulsar, was initially detected through its 'gia nt' pulses, 
not through its t ime-averaged pulsed emission dStaclin &| 
Reifcnstcin 1968] ). Giant pulses are individual pulses with 



amplitudes much grea ter than the mean pulse amplitude 
( Lundgrcn ct al. 1995 ). As we might expect other young 
pulsars to behave similarly to the Crab, we also searched 
for aperiodic, dispersed pulses from the XPS. We recorded 
all pulses with amplitudes above a signal-to-noise thresh- 
old of 4 for each dedispersed time series, enchancing our 
sensitivity to broadened pulses by repeating this thresh- 
olding with different levels of time-series smoothing. Given 
our 4ct threshold, we were sensitive to all single pulses with 
flux densities greater than 25/ \/w m s Jy where w m s is the 
pulse width in milliseconds, at 327 MHz and 1.0/y/w ms 
Jy at 1435 MHz. Our smoothing algorithm optimized the 
search for pulse widths up to 15 ms. (For comparison, 
the widths of the Crab pulsar's giant pulses range from 
roughly 0.2 to 0.4 ms.) We found no evidence for isolated 
pulses from the XPS. The result of this analysis, along 
with the result for PSR B0355+54, is shown in Figure 2. 
For comparison, the Crab pulsar emits roughly 80 pulses 
per minute with flux densities greater than 100 Jy at 327 
MHz and 8 pulses per minute with flux densities greater 
than 50 Jy at 1435 MHz. 
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Fig. 2. — The number of isolated, dispersed pulses at 327 MHz 
above a 4<r threshold is plotted vs. DM for a 5-minute observation 
of PSR B0355+54 (left) and a 50-minute observation of the XPS 
(right). While a peak at the 57 pc cm' 3 DM of PSR B0355+54 is 
obvious, there is no evidence for such an excess from the XPS. The 
peak at low DM in the right plot is consistent with interference, and 
is similar to the number of counts observed in the consecutive blank 
sky pointing. I n the absence of radio frequency interference and/or 
a pulsar, we would expect a flat distribution. 



4. RESULTS 

To calculate an upper limit for pulsed radio emission 
from the XPS, we must account for the excess system tem- 
perature due to the em ission from C as A itself. The Cas A 
remnant is 5' in size (Green 2000). Hence, the ratios of 
our beam area to remnant area are 0.003 and 0.0002 at 
327 and 1435 MHz. Given Cas A 's measured flux density 
and spectral index ( Green 2000D , we calculate remnant 
flux densities of 6432 and 2060 Janskys at 327 and 1435 
MHz, contributing to an increase in system temperature 
of a modest 40 K at 327 MHz and a mere 1 K at 1435 MHz 
for the phased array. However, while the phased array re- 
solves out emission from the nebula, each single dish will 
see the whole remnant, leading to an increase in system 
temperature of 475 K at 327 MHz and 230 K at 1435 MHz. 
By comparison, a single dish with the same collecting area 
as the VLA would suffer increases of system temperature 
of 13000 and 6200 K at 327 and 1435 MHz. 

We calculate the minimum detectable flux density o f our 
search using the expression ( pordes fc Chcrnoff 1997] ) 



iS,, 



G^/N pol AuAtN FFT J2fL\ R, 



(1) 



where rj is the signal-to-noise threshold used, T sys is the 
system temperature, G is the telescope gain, N po i is the 
number of polarization channels, Av is the total band- 
width, At is the sample interval, TVfft is the number of 
points in the FFT, is the optimal number of harmon- 
ics with which a pulsar can be detected, and Ri < 1 is the 
ratio of the Zth harmonic to the fundamental. Nh will de- 
pend on the pulsar's duty cycle, or ratio of pulse width to 
period, with narrower pulses producing more harmonics. 
For this calculation, we take the intrinsic duty cycle to be 
the minimum of 0.3 and 0.03(P ms /1000.) ~ 0,5 , where P ms is 



the pulse period in milliseconds (Biggs 1990). Pulses may 
then be additionally broadened by propagation effects, in- 
cluding scattering and dispersion. The Taylor & Cordes 
(1993) model predicts 1.7 ms and 3 /xs of pulse broadening 
due to scattering in the direction of the XPS and for a 
distance of 3.4 kpc at 327 and 1435 MHz, respectively. 



Figure 3 shows the dependence of our search sensitiv- 
ity on DM and period at the two observing frequencies, 
accounting for the effects of scattering and dispersion on 
the intrinsic duty cycle of the pulsar. Assuming values for 
period and DM of 30 ms and 75 pc cm~ 3 , we calculate 
7a upper limits of 30 and 1.3 mJy at 327 and 1435 MHz. 
As shown in Figure 3, these limits will decrease for longer 
period, smaller duty cycle, and lower DM. There are no 
hardware limitations to our maximum detectable period, 
determined solely by the length of the FFT. 

While there is no infrared or optical evidence f or a bi- 
nary companion (van den Bcrgh & Pritchct 1986), accel- 
eration effects from such a companion would cause pulse 
smearing, decreasing our sensitivity. We estimate the 
signal-to-noise degradation introduced for various orbital 
systems by taking a Gaussian pulse and applying orbital 
smearing by iteratively calculating emission times across 
the orbit. Smearing times due to time constant, DM, 
and scattering are calculated and applied to the orbital- 
smeared pulse. We calculate the best harmonic sum for 
the resultant waveform and compare this to the best har- 
monic sum for a waveform with no orbital smearing to 
calculate the loss in signal-to-noise. For a 30 ms period 
and duty cycle of 0.2, a 1 solar mass companion in an 8 
hour circular orbit with an inclination angle (i.e. angle 
between the orbital plane and our line of sight) of 30°, our 
sensitivity decreases by roughly a factor of 2, averaging 
over all orbital phases, for a 16 Mpoint FFT. In the case 
of our 8 and 4 Mpoint FFTs, these same binary param- 
eters would decrease our search sensitivity by factors of 
2.4 and 2.8 over the original 16 Mpoint sensitivity to an 
unaccelerated pulsar. Although the nominal sensitivities 
are less for these shorter FFT lengths, the shorter FFTs 
would allow us to detect pulsars in fast binaries where the 
binary smearing during a longer FFT is comparable to or 
greater than the pulse period. The sensitivity would be 
further decreased for smaller spin periods, faster orbits, 
and more massive companions. If the orbit was eccentric, 
the net smearing, averaged over all orbital phases, would 
be less, and the resultant sensitivity would be higher. Ac- 
cretion from a binary companion could completely quench 
the radio emission. 
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Fig. 3. — Contours of pulsed search 7<j sensitivity vs. P and DM 
for searches at both 327 and 1435 MHz. The thickest (rightmost) 
line corresponds to a sensitivity of 16 mJy at 327 MHz and of 0.8 
mJy at 1435 MHz. Minimum detectable flux densities increase by 
a factor of 2 with eac h subsequently thinner line. The contours are 
not smooth because of the finite spacing of our trial DMs. 

Interstellar scintillation modulations of the flux of 
the XPS are mostly negligible during our observations. 
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The predicted scintillation bandwidth and timescale (see 
Cordes & Rickett 1998), assuming a uniform medium and 
a velocity of 50 km s _1 are approximately 0.1 kHz and 
30 seconds at 327 MHz and 0.6 MHz and 150 seconds at 
1435 MHz, in both cases too small to significantly affect 
our quoted upper limits. If this velocity is higher than 50 
km s -1 , as a you ng pulsar's velocity is likely to be ( Cordcij 
& Ch ernoff 1998| ), the predicted scintillation timescale will 
be smaller and will therefore be even less of a factor in our 
analysis^ 



We have also used VLA images of Cas A to search for a 
radio continuum point-source at the position of the XPS. 
At 1.3 GHz and 4.4 GHz, we find 5<r upper limits on the 
flux density from the XPS of 40 mJy and 6 mJy respec- 
tively. While not as constraining as those derived from our 
VLA and Arecibo pulsed searches described above, these 
continuum upper limits are completely independent of the 
period, dispersion and scattering of any radio pulsations 
from this source. 



5. CONCLUSIONS 

We have calculated upper limits of 30 and 1.3 mJy at 
327 and 1435 MHz for the phase-averaged pulsed flux den- 
sity from the XPS in Cas A. We did not detect any single, 
dispersed pulses above our sensitivity of 25 and 1.0 Jy (as- 
suming a pulse width of 1 ms) at 327 and 1435 MHz. We 
now compare our upper limit to those from previous Cas A 
searches. Davies & Large (1970) carried out an unsuccess- 
ful 408-MHz search for single pulses from Cas A. Seiradakis 
& Graham (1980) quoted a 1420-MHz upper limit of 1.6 
mJy for pulsed emission from Cas A. However, their limit 
does not account for the additional system temperature 
contributed by the nebula, significant for their beamwidth 
of 9'. More recently, Woan & Duffett-Smith (1993) pub- 
lished an upper limit of 80 mJy at 408 MHz and Lorimer et 
al. (1998) published an upper limit of 46 mJy at 606 MHz. 
For a typical pulsar spectrum, our 1435-MHz upper limit 
is an order of magnitude better than these previously pub- 
lished limits, due to the known position of the XPS and 
the power of the VLA to resolve out the excess emission 
from the nebula. 

Our upper limit of 1.3 mJy at 1435 MHz translates to a 
luminosity of 15 mJy kpc 2 for a neutron star at a distance 
of 3.4 kpc. This is lower than the mean 1400-MHz lumi- 
nosity of 30 mJ y kpc 2 estimated for young pulsars ( Fraij 
& Moffctt 1993) and is lower than the luminosity of any 
known radio pulsar with characteristic age less than 10 4 
years (see Figures 1 and 2 of Brazier & Johnston 1999). 
This suggests that the XPS in Cas A is either not a pul- 
sar, is a radio-quiet pulsar, or is a radio pulsar beamed 
away from us. While evidence suggests that pulsar beams 
are wider at short periods (see e.g. Biggs (1990)) and that 
they are therefore likely to intersect the Earth's line of 
sight, individual pulsars show mixed properties. For in- 
stance, the Crab pulsar shows narrow radio components 
that may be part of a large beam, but the profile of another 
young pulsar in th e Large Magellanic Clou d, PSR B0540- 
69, is very broad ( Manchester et al. 1993| ). Furthermore, 
it is difficult to disentangle true beam width from viewing 
angle for these pulsars. It is certainly possible that, like 
the Geminga pulsar, the XPS shows no radio pulses but 
is an X-ray, and perhaps gamma-ray, pulsar. Although 
there is no EGRET point source near the position of the 



XPS, gamma-ray searches with future, more sensitive in- 
struments, could be successful. Further measurements of 
the spectrum and time variability of the XPS will also help 
determine which of the above scenarios is the case. 

This search adds to the growing body of evidence that 
there arc manifestations of young neutron stars other 
than the standard Crab- like radio pulsar. Although it 
is widely accepted that radio pulsars are born in super- 
novae, fewer than 10% of searches targeted at sup ernova 
remnants h a ve found young (age < 25 kyr) pulsars (Kasp: 
et a. 1996|, |Gorham et al. 1996;, |Lorimcr et al. 1998). 



Some of this may be due to scattering at low frequencies 
in the Galactic plane and the decreased sensitivity due 
to the bright remnants, but it also hints at the unrec- 
ognized diversity of the neutron star population. While 
some neutron stars are detected as radio pulsars, others 
may be born as anomalous X-ray pulsars, soft gamma-ray 
repeaters, quiescent neut ron stars, or other classes of ob- 
jects not yet discovered (Kaspi 200C, |Gotthelf fc Vasisht 



2000 ) . The pathways to and percentages of these different 



classes of objects are not yet clear. However, future Chan- 
dra detections of compact X-ray sources, and large-scale 
radio pulsar surv eys, like the Parkes Multibeam Survey 
(Lyne et al. 2000), will help to answer these questions. 
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